Reconstruction of the second layer of Ag on Pt(lll) 
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The reconstruction of an Ag monolayer on Ag/Pt(lll) is analysed theoretically, employing a 
vertically extended Frenkel-Kontorova model whose parameters are derived from density functional 
theory. Energy minimization is carried out using simulated quantum annealing techniques. Our 
results are compatible with the STM experiments, where a striped pattern is initially found which 
transforms into a triangular reconstruction upon annealing. In our model we recognize the first 
structure as a metastable state, while the second one is the true energy minimum. 
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I. INTRODUCTION 

One of the most important physical processes that 
must be properly understood in order to manufacture 
the nanoscale devices is surface reconstruction, i.e., the 
morphology changes in the epilayer of a material as a 
result of reduced coordination, either in homoepitaxial 
or heteroepitaxial systems. The lattice parameter of the 
epilayer is usually different from that of the bulk system, 
giving rise to surface stress. This stress can be relieved in 
a variety of ways, inducing the appearance of interesting 
patternsi^. 

In this work we analyse the reconstruction of the sec- 
ond monolayer (ML) of Ag deposited on a clean surface 
of Pt(lll). The experiments of Brune and coworkers^ 
show that, while the first ML is pscudomorphic (i.e., 
copies the structure of Pt(lll)), the second monolayer 
reconstructs in a striped pattern, which upon annealing 
at 800 K, gives rise to an intrincate trigonal network, 
reminiscent of a kagome lattice. Along (111) direction, 
silver presents a lattice parameter of aA g = 2.95 A, which 
is 4% higher than that of platinum, apt = 2.83 A. 

Henceforth, we will refer to the system composed of 
a complete first monolayer (ML) of Ag on Pt(lll) as 
Ag/Pt(lll). Ag atoms in the second ML on that struc- 
ture can stay in two types of local minima, both arranged 
in triangular networks, which are denoted as FCC and 
HCP sites (Figure [T]). Earlier work&£ has shown that 
the two types of sites have very similar energies, being 
the FCC site favoured by 3 meV. In the middle point of 
the segment connecting a FCC and a HCP site there is 
a bridge saddle-point site. The energy of an Ag atom in 
the second ML at the bridge site constitutes the energy 
barrier to jump from one type of minimum to the other. 
It was estimated by a theoretical fitting of experimental 
data on island densities^ to be around 60 meV. As an 
important final remark, the adatom-adatom interaction 



was shown to decay slowly, with a repulsive ring past 
the short-range which may play an important role in the 
kinetics^. 

In the present work we model the second monolayer 
of Ag on Ag/Pt(lll), referred to as Ag/Ag/Pt(lll). 
We use a generalization of the renowned 2D Frenkel- 
Kontorova (FK) model^. This model has provided ac- 
curate and insightful analysis of similar reconstruction 
processes, such as Pt(lll)^ or the first layer of Ag on 
Pt(lll)ii. Our model differs from the standard FK in 
the following aspects. First of all, the substrate potential 
is vertically extended, i.e., it becomes a function of x, y 
and z. This way, the film is allowed to relax the stress 
by (small) vertical displacements of the atoms in the top- 
most Ag layer—. A second difference is that the number 
of neighbours of any atom is not fixed beforehand, i.e., 
no topology is assumed. Thus, all kinds of lattice defects 
are allowed in the surface structure. 

The substrate potential for the extended FK model is 
obtained by choosing a physically sensible general form 
and obtaining its parameters using ab initio density func- 
tional theory (DFT). For the film potential, we use a 
Morse type pair potential. The obtention of the min- 
imum energy configurations of the extended FK model 
is a complex problem, which has been approached using 
the recently developed replica-pinned quantum annealing 
(RPQA)^. 

The minimum energy configurations of the extended 
FK model are assumed to be (local) equilibrium atomic 
patterns of the reconstructed surface. Of course, this ap- 
proach neglects the kinetic effects which may be relevant 
in the determination of the experimental configurations. 
We will provide some comments about these effects in 
section IIVI 

This article is organized as follows. SectionlHl discusses 
the DFT calculations and the obtention of the extended 
FK model. The results of the numerical minimization 
are given in section Hill while its physical interpretation 
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FIG. 1: Geometry of the Ag/Pt(lll) system: top sites are 
shown as the vertices of the triangular network; FCC, HCP 
and bridge sites are shown as empty circles, empty squares 
and filled dots respectively. A unit cell is marked with a 
thicker line. 



is provided in section HV1 along with the most promising 
lines for further work. 



II. MODEL 

The behavior of an Ag monolayer on the Ag/Pt(lll) 
system is simulated by an extended FK model, defined 
by the following functional form 

E({n}) = E + E VKk< " r iD (!) 

where rj = (xi, yi,zi) is the 3D position vector of the i th 
Ag atom, V s denotes the substrate potential, Vf is the 
film potential, and the sum over pairs (i,j) is extended 
over neighbours within a certain cutoff distance. This 
cutoff distance is determined so as to include the ring of 
first NN atoms. The vertical extension of the substrate 
potential^ allows for stress relaxation via small vertical 
displacements along z-direction. 

The appropriate intensive magnitude to minimize is, of 
course, the energy per unit area. The unit cell is shown 
by a bold line in Figure |T]). The atomic density at the 
surface is, therefore, another relevant parameter to be 
obtained in our method. 

Let us review concisely the geometry of our model. 
There are four types of distinguished points on the sub- 
strate lattice, which are known as FCC, HCP, top and 
bridge, as shown in Figure [TJ The local equilibrium posi- 
tions for an atom are the FCC and HCP, being the first 
slightly favoured energetically. Bridge sites are saddle- 
points (usually denoted with the letter p), and top po- 
sitions, so called because they stand exactly above the 
substrate atoms, are strongly disfavoured energetically. 

A. Density Functional Theory calculations 

The parameters in the substrate potential arc obtained 
using ah initio density functional theory (DFT) calcula- 



tions with the pscudopotential approach and plane waves 
basis, as implemented in Quantum-Espresso^. We use 
generalized gradient approximation (GGA) for the ex- 
change correlation interaction with the functional pro- 
posed by Perdew, Burke and Ernzerhofi^. A kinetic 
energy cut-off of 367.2 eV and a higher cut-off of 2940 
cV is used for the augmentation charges introduced by 
the ultrasoft pscudopotcntiali^. To improve the conver- 
gence, a gaussian smearing of width 0.68 eV was adopted. 
Brillouin zone integrations for the (1 x 1) surface cell, 
shown by a thick line in Figure [TJ were carried out using 
a (15 x 15 x 1) mesh of fc-points. 

We obtained the bulk lattice parameter for Ag and 
Pt as 4.17 and 4.0 A respectively, which give a misfit of 
4.25%. These values compare well with the experimetal 
values, 4.09 and 3.92, with a misfit of 4.33%, and the 
previous theoretical calculations of 4.19 A and 4.01 A^. 
According to our values, the nearest neighbour (NN) dis- 
tance for Ag and Pt on the (111) surface are 2.95 and 
2.83 A respectively. 

B. Substrate potential 

The substrate potential V s (r) must share the periodic- 
ity of the substrate lattice, i.e., it will take the form of a 
Fourier series^, 

V s (r,z) = V + Y,V G (zy G - r , (2) 

G 

where G are 2D- vectors of the reciprocal substrate lattice 
and r = (x, y). This series is truncated arbitrarily to con- 
tain the first two shells of G vectors, of length 47r/v / 3apt 
and 47r/apt, separated by 27r/3 rad, with apt = 2.83 
A. The dependence on z is restricted to the coefficients 
Vg{z). When we expand equation (J5J in the Vfe's using 
the first two shells of G vectors and apply the symmetries 
of the surface, we get an equation with four unknowns^. 
In order to obtain these four unknowns, we need to know 
the values of the substrate potential (V s ) at four non- 
equivalent positions on the surface, which we choose to 
be the FCC, HCP, bridge and top sites. Thus, we express 
the substrate potential V^(r) in terms of Vfcc, Vhcp, 
V p , V t0 p- These values are calculated within DFT for a 
monolayer of Ag on a Ag/Pt(lll) slab. The Ag/Pt(lll) 
slab is simulated using a symmetric slab of seven Pt lay- 
ers and one Ag layer adsorbed on both sides of Pt. We 
first calculate the equilibrium heights by relaxing the Ag 
monolayer on Ag/Pt(lll) along the z-direction at the 
four sites. We use a vacuum width of 8 layers (18.47 A) 
between the two Ag/Ag/Pt(lll) slabs with the in-plane 
lattice parameter of 2.83 A. Since the Vg's are functions 
of z, we calculate the energy of the second ML of Ag at 
various values of z by keeping the monolayer at different 
heights from the substrate. In this way, four functions of 
z are obtained: i>k(z), with k G {FCC, HCP, p, top}, by 
fitting the DFT results to a Morse-like form: 
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FCC 


HCP 


bridge 


top 


Zeq (A) 


2.45926 


2.45977 


2.54478 


2.80834 


E min (meV) 





3.638 


35.2784 


174.066 



TABLE I: Equilibrium heights (z eq ) and corresponding ener- 
gies to create an stacking fault in the second ML of Ag with 
respect to lowest energy FCC site, as obtained from DFT. 



*(A) 


E (FCC) 


E (HCP) 


E (bridge) 


E (top) 


2.27447 


21.3928 


26.7308 


69.2716 


326.563 


2.64406 


17.4284 


19.2576 


42.4184 


192.134 


2.73646 


34.4624 


36.7064 


57.6368 


178.520 


3.34878 


252.375 


253.416 


254.735 


269.606 



TABLE II: Energy (in meV) of an Ag atom on each of the four 
types of sites, at some values of z, with respect to the energy 
of the atom at its equilibrium height on FCC, as obtained 
with DFT. 



Ufc(z) 



A 



"*~ V k ,o 
(3) 

where = /ifc/2. The substrate potential is now well 
defined from ab initio calculations. 

If we evaluate the substrate potential at the equilib- 
rium height of the FCC position, we should obtain the 
chemical potential for the Ag/Ag/Pt(lll) system, fi, i.e., 
the energy required to take a single Ag adatom from its 
equilibrium position on the surface to infinity. We have 
estimated \i from DFT calculations as three times the 
binding energy of two Ag atoms in bulk. Thus, fi = 1.251 
cV. 

Table U gives the equilibrium heights calculated from 
DFT at each type of site, along with the energy required 
to create a stacking fault relative to that of the FCC po- 
sition. The value for the energy difference between the 
FCC and the HCP sites is consistent with the previous 
values given in the literature, about 3 meV. The energy 
barrier between them is lower in our case, about 35 meV, 
as compared to the 60 meV of other sources^. The rea- 
son for this discrepancy is that they have done calcu- 
lations on a single adatom adsorbed on FCC or bridge 
sites, whereas we use one full monolayer of Ag on Ag/Pt. 
Using a single adatom on FCC or bridge sites would ob- 
viously increase the binding energy of the adatom due to 
its lower coordination on the surface with respect to a 
full monolayer. 

Some DFT data for the dependence on z of the sub- 
strate potential are given in Table HH A plot of V s (r) — fj, 
is shown in Figure [2j We would like to remark that the 
fit of V s (r) includes more points than those shown in ta- 
ble [IT] (between 5 and 8 per site type), but only those are 
shown because they are computed for equal values of z. 



y(A) 




z (angstroms) 

FIG. 2: (Color online) The substrate potential V s (r) — fi. Up- 
per: contour plot of V s (x, y, zf cc ) — fi calculated at the equi- 
librium height of the FCC sites. Lower: vertical dependence 
of the substrate potential, for the FCC, HCP, top and bridge 
sites. 



C. Film potential 

The film potential Vf (r) (also known as surface poten- 
tial) simulates the interaction among Ag atoms in the 
Ag layer on Ag/Pt(lll) surface. We assume it to be an 
isotropic pair potential of Morse type. 



V f (r) = a 



-a t (r-b) 



(4) 



In this equation, b is the equilibrium distance, oq the 
depth of the potential well and ai the curvature around 
the minimum. However, Vf(r) is rather difficult to obtain 
from DFT calculations due to the heteroepitaxial nature 
of the problem. The closest approximation would be to 
calculate it for Ag atoms on a Ag(lll) surface. However, 
the nearest neighbour (NN) distance in the Ag layer on 
Ag/Pt(lll) surface is 2.83A, while the NN atomic dis- 
tance in the Ag layer on Ag(lll) surface is 2.95A. There- 
fore, as one would expect, Ag layer on Ag/Pt(lll) sur- 
face is under compressive stress of 0.198 eV/A 2 and Ag 
layer on Ag(lll) surface is under tensile stress of 0.091 
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eV/A 2 . Due to this difficulty, we have fixed the value 
of b to 2.95A, i.e., the lattice parameter of Ag, and run 
simulations for different values of do and a±, as we will 
see in the next section. 

In order to calculate the total film energy, we only con- 
sider pairs of atoms closer than a certain cutoff distance 
^cutoff =6A. As a technical side remark, we added a tiny 
constant to the potential in order to avoid a discontinuity 
at that distance. Therefore, the topology of the epilayer 
is not fixed beforehand: the number of neighbours of a 
given atom is not restricted. 

The Ag-Ag pair potential has a complex structure at 
longer distances, which we have neglected in this work. 
Other authors^ have shown the effect of Friedel-like oscil- 
lations at long distances, recognizing a ring of repulsion 
in the pair potential. This ring is beyond our cutoff dis- 
tance and we will assume it to be negligible in order to de- 
termine the equilibrium configuration. Its effects on the 
kinetics of the system might be much more pronounced, 
and will be discussed in section llVl 
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FIG. 3: Top view of the striped reconstruction of the 
Ag/Ag/Pt(lll) surface, on a 25 x 4 unit cell. Notice the 
wavy pattern, with a wavelength of « 25ap t . Heights are de- 
noted by the shades of gray: the darker the atom, the lower 
its height. 

Detailed results and comparisons will be provided else- 
where. 

Of course, it is impossible in these problems to be com- 
pletely sure that the best minimum obtained is the true 
global minimum, but after an exhaustive search and co- 
incidence of diverse methods a large degree of confidence 
can be assumed. 



B. Results 



III. OPTIMUM CONFIGURATIONS 
A. Numerical optimization 

Once the extended FK Hamiltonian has been deter- 
mined, i.e., the substrate and film potential have been 
found, the main task is to obtain the positions of the Ag 
atoms which minimize the total energy, i.e., its equilib- 
rium configuration. 

Global optimization of a potential energy surface 
(PES) with a large amount of local minima is a highly 
non-trivial task. Simulated thermal annealing is one of 
the most popular non-biased general-purpose optimiza- 
tion methods, in which thermal fluctuations allow the 
system to escape metastable states in order to achieve the 
real global minimum. A recent innovation is that of simu- 
lated quantum annealing, in which quantum fluctuations 
collaborate with thermal fluctuations in this tas k 18 ' 19 . 
Using path integral concepts, the original system is re- 
placed by a set of replicas exploring the same PES, linked 
among themselves with springs of zero natural length, 
at a fictitious finite temperature. Replica-pinned quan- 
tum annealing (RPQA) is a technical variant which has 
provided the best non-biased results to date for the de- 
termination of the global minimum energy structures of 
Lenard-Jones clusters. The basic idea is to update all 
replicas excluding the one which has reached the best 
energy so far. The details can be found in the original 
article by Gregor and Car—. 

The application of the RPQA to the obtention of the 
global minimum of FK models is one of the main tech- 
nical improvements developed in this work. Both in ID 
and in 2D, the efficiency of the method is higher than 
the case of simulated thermal annealing. RPQA tends to 
explore more thoroughly the area around the best replica 
so far, thus removing defects faster than other methods. 



An exhaustive search for the minimum energy configu- 
rations of the extended FK model, as defined in equation 
[U has been carried out using the procedure described in 
the previous section. For all the simulations, we use peri- 
odic boundary conditions (PBC), with different sizes for 
the system cell. A few important parameters were sys- 
tematically varied. On one hand, the parameters ao and 
a±, i.e., the depth and curvature of the film potential, 
for which ab initio methods have not provided a reliable 
estimate. On the other hand, the atomic density of the 
surface, i.e., the number of particles per unit area. 

Among the huge number of local minima, we have se- 
lected three of them because of their special physical rel- 
evance. They are local minima of the PES for all reason- 
able values of a (l and a\. 

(a) Unreconstructed configuration. The homoge- 
neous configuration in which all Ag atoms stay in 
the FCC positions on Ag/Pt(lll). The surface is 
under compressive stress. 

(b) Striped configuration. In this configuration the 
stress is relaxed by placing some particles in HCP 
positions, as shown in Figure ([3]) (25 x 4 unit cell 
with 24 x 4 Ag atoms). The increase of substrate 
energy is compensated by the decrease in the film 
energy due of the stretching of the horizontal links. 
The wavy pattern established has a preferred wave- 
length, of approximately 25ap t . Thus, there is a de- 
crease of density as compared to the first Ag ML, 
of 4%. 

(c) Triangular configuration. Stress is relaxed in 
an isotropic way, creating lines of defects following 
the symmetry axes of the substrate, as it is shown 
in Figure 3] The density of this configuration de- 
creases by 7.8% with respect to the unreconstructed 
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FIG. 5: Phase diagram which shows the global minimum of 
the PES as a function of the parameters of the film potential, 
an and oi. 
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FIG. 4: Top view of the triangular reconstruction of 
Ag/Ag/Pt(lll) surface. Heights are denoted by the shades 
of gray: the darker the atom, the lower its height. 



case. In order to allow the relaxation of the surface 
bonds, some atoms must stay close to top sites. 

Which of them is the global minimum of the PES de- 
pends on the values of ao and oi. Figurc[5]dcpicts a phase 
diagram in which the stable configuration is shown as a 
function of these two parameters. 

The striped configuration is reported experimentally 
by Brune and coworkers^. According to their work, un- 
der annealing at 800 K, it relaxes to another configu- 
ration which is similar to our triangular configuration. 
It is well known that, under annealing, the system is 
much more likely to find the global energy minimum of 
the PES. Therefore, the experimental data imply that 
the triangular structure is, indeed, that global minimum, 
while the striped configuration is a long-lived metastable 
state, i.e.: a local minimum. The reasons for its long 
lifetime are a large energy barrier and a small difference 
in energies between both configurations. 

Therefore, the numerical data are compatible with the 
experimental results in as much as the values of ao and 
a\ are close to the curve which separates the striped and 
the triangular regions in the phase diagram of Figure 
([5|). The precise point within that curve can not be de- 
termined from our calculations. 



IV. DISCUSSION AND FURTHER WORK 

Our theoretical model shows from ab initio calcula- 
tions, that the Ag/Ag/Pt(lll) can reconstruct into a 
striped or a triangular configuration for some choices of 
its free parameter (ao and a\ in the film potential). Both 
the equilibrium and the metastability situations can be 
explained satisfactorily. 

The fact that the striped configuration is preferred be- 
fore annealing requires further clarification. It must have 
a larger energy per unit cell than the triangular struc- 
ture, but kinetic effects may make it more likely to ap- 
pear. A full theoretical explanation of this fact remains 
as an open task, but we may already make some remarks 
with the available information. At room temperature, Ag 
adatoms on the Ag/Pt(lll) surface will not distinguish 
clearly between HCP and FCC sites. The repulsion ring 
described by Fichthorn and Schcffier— explains the forma- 
tion, at low coverages, of a high amount of small islands 
surrounded by an exclusion area which makes it difficult 
for other adatoms to join. With the energy difference of 
3 meV between FCC and HCP sites, we can predict that, 
at room temperature, roughly 50% of these islands will be 
on each type of site. Islands will merge slowly, and some 
domain boundaries will appear. At moderate coverages, 
the most likely configuration for these domain bound- 
aries is a set of parallel lines, giving rise to the striped 
configuration. The complex network which is character- 
istic of the triangular configuration will be more unlikely, 
although it leads to a lower energy for a complete layer. 

There are some differences between the global mini- 
mum with triangular symmetry obtained in this work 
and the configuration reported by Brune and coworkers. 
In the later structure, no atoms sit on bridge or top sites. 
With our PES, the energy per unit cell of that configu- 
ration is higher than the energy of the unreconstructed 
surface. When we use this structure as the starting con- 
figuration for our quantum annealing minimization al- 
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gorithm, the system relaxes to our triangular structure. 
A thorough search has been performed in the parame- 
ter space in order to obtain the structure of Brune and 
coworkers as the global minimum of the PES, without 
success. This could mean that, we will have to abondon 
the Frenkel-Kontorova approximation in order to obtain 
the experimental configuration as a global minimum of a 
PES. The assumption which is most likely to fail is the 
independence of the film potential from the coordination 
number. The lower coordination of some sites increases 
the strength with which they bind to their neighbours. 
This effect is taken into account in other more complex 
models, such as the glue models. 

Other systems which are of interest, where our ver- 
tically extended Frenkel-Kontorova model can provide 
insight are the reconstruction of similar metallic sur- 
faces, such as Cu/Ru(llll). Surface growth in semi- 
conductors, such as the Stranski-Krastanov growth mode 
in InAs/GaAs will require major modifications, since 
isotropy of the film potential might not be an acceptable 
approximation. On the other hand, the vertical exten- 



sion could prove much more useful in this case, where 
the film tends to really curve. 

On the computational side, this work applies quan- 
tum annealing methods to a physical problem with 
experimental comparison, i.e., not just a benchmark 
problc m 13 i 18 i 19 . Of course, quantum annealing should 
not be considered as the panacea for all optimization 
problems. It is just another tool in our numerical tool- 
box, which can prove useful in a variety of cases. 
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